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ABSTRACT 

The radial entropy profile of the hot gas in clusters of galaxies tends to follow a 
power law in radius outside of the cluster core. Here we present a simple formula 
giving both the normalization and slope for the power-law entropy profiles of clusters 
that form in the absence of non-gravitational processes such as radiative cooling and 
subsequent feedback. It is based on seventy-one clusters drawn from four separate 
cosmological simulations, two using smoothed-particle hydrodynamics (SPH) and two 
using adaptive-mesh refinement (AMR), and can be used as a baseline for assessing 
the impact of non-gravitational processes on the intracluster medium outside of clus- 
ter cores. All the simulations produce clusters with self-similar structure in which the 
normalization of the entropy profile scales linearly with cluster temperature, and these 
profiles are in excellent agreement outside of 0.2r2oo- Because the observed entropy 
profiles of clusters do not scale linearly with temperature, our models confirm that 
non-gravitational processes are necessary to break the self-similarity seen in the simu- 
lations. However, the core entropy levels found by the two codes used here significantly 
differ, with the AMR code producing nearly twice as much entropy at the centre of a 
cluster. 
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1 INTRODUCTION 

Purely gravitational structure formation ought to produce 
clusters of galaxies with nearly self-similar structure, whose 
X-ray luminosity Lx scales with temperature Tx as Lx oc T x 



jKaiserl Il986f) . Clusters created in hydrodynamical simu- 
lations with cosmological i nitial conditions indeed f ollow 
this scaling relation (e.g., iNavarro et alJ Il995l : lEke et ail 
1998), but observed clusters do not. Instead they follow 



a relation closer to Lx oc T 2 ' 8 (Edge & Steward 


1991 


David et alJll995l: 


Markevitch 1998; Allen & Fabian 


1998 


Arnaud & Evrardl 


1999). Somehow the non-gravitational 



cooling and heating processes associated with galaxy for- 
mation intervene to break the expected self-similarity, with 
consequences that are more severe in lo w-temperature 
clusters than in high -temperature clusters jKaiserl Il99ft 
lEvrard fc Henrvlll99ll . see Voit 2005 for a recent review). 

The scaling behavior of the Lx-Tx relation can largely 
be understood in terms of radiative cooling and the feed- 
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back it triggers tVbit fc Brvanl l200ll) . Gas that can cool 
within a Hubble time must either condense, forming stars 
or cold baryonic clouds, or it must be reheated somehow, 
probably by supernovae or AGN activity triggered by the 
condensing gas. Because ~ 85% of the baryons associated 
with a massive cluster appear to be in the hot intraclus- 
ter medium (ICM), a certain amount of supernova or AGN 
feedback seems necessary to prevent too many of the baryons 
fro m condensing during t he formation of a cluster's galax- 
ies JWhite fc Frenklll99ll : rBalogh et af]l200ll) . However, the 
Lx-Tx relation itself is relatively insensitive to the total 
amount of feedback energy, as long as the energy input into 
the reheated gas is sufficient to keep it from cooling again 
jBoreani et alJl200l iKav et alJliooa: iTornatore et alJl2003l: 
IValdarninill2003l) . Thus, one needs to look beyond the Lx- 
Tx relation in order to assess the thermodynamic impact of 
supernovae and AGNs on the state of the ICM. 



One good place to look for more information is in the 
spatially-resolved entropy profiles of clusters and groups, 
which preserve a record of the cooling and heatin g processes 
responsible for similar i ty breaking in clusters jVo it et alJ 
120021: IVoit et alJ I200&1: iKavl 12004: IKav et alJ 12004) . Spe- 
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cific entropy, represented in this paper by the quantity 
K = Trie 2//3 , where n e is the electron density, is more closely 
tied to the thermodynamic history of a cluster than is tem- 
perature, because the thermal energy of heated gas can be 
converted into gravitational potential energy as the heated 
gas expands in the confining potential well. 1 Under certain 
circumstances, a large amount of energy input produces only 
a small rise in th e luminosity- weighted temperature Tx (e.g., 
IVoit et al] 12002ft . The specific entropy of the ICM, on the 
other hand, always rises when heat energy is introduced and 
always falls when radiative cooling carries heat energy away. 

The entropy profiles of clusters and groups can now be 
measured out to a significant fraction of the scale radius T2oo, 
within which the mean mass density is 200 times the crit- 
ical density. Those measurements show that entropy levels 
in the cores of clusters, where the Lx-Tx re lation is deter- 
mine d, scale as Ko.i = if(0.1r2oo) oc T^/ 3 iPonman et alJ 
2003), as expected if radiative cooli ng and associated feed - 
back govern the core entropy level J Voit fc Ponmanll2003ft . 
More surprisingly, entropy measurements at larger cluster 
radii are hinting that this scaling relation applies to the en- 
tire entropy profile. Deep XMM-Newton observations of five 
clusters whose temperatures span a range of ~3.5 show that 
the scaled pr ofile T^ 2 ^ 3 K(r /r2oo) is independent of cluste r 
temperature JPratt fc Arnaudl2003l : lPratt fc Arnaudl20o3) . 
Likewise, an analysis of lower-quality data on a larger num- 
ber of clusters also suggests that K(r/r2oo) oc at the 
scale radius r^oo, within whic h the mean matter d ensity is 
500 times the critical density jPonman et al]l2003ft . 

Rather than totally breaking the self-similarity of clus- 
ters, galaxy formation appears to alter the power-law scal- 
ing of K(r/r2oo) with Tx without appreciably changing the 
overall shape of the entropy profile. Exactly how heating and 
cooling would conspire to produce such a shape-preserving 
shift in the normalization of an intracluster entropy profile 
is unknown. One possibility involves smoothing of the inter- 
galactic medium by supernovae or AGN energy input prior 
to accretion, which lowers the mass-weighted mean density 
Pace of the infalling gas. Because the amount of entropy gen- 
erated in that gas when it passes through accretion shocks 
of velocity v aC c is K^ cc ~ t'acc/Pacc, smoothing of gas that 
would otherwise be bound to accreting su bhalos boosts the 
post-accretion entropy level of the ICM llVoit et al ] Li M it 
iPonman etld1l2003l: iBorgani et al.ll2005|). An other possibil- 
ity, illustrated in simulations by Kav (2004), involves un- 
compensated cooling, which allows high-entropy gas to sink 
to smaller radii as the core gas condenses, but pure cooling 
does not appear to reproduce the K oc T^/ 3 relation at rsoo- 

A proper analysis of the observations to determine the 
true source of the entropy boost requires knowing what the 
baseline entropy profile of a cluster would be like in the ab- 
sence of galaxy formation. To that end, this paper compares 
the results of four different hydrodynamical simulations of 
purely gravitational cluster formation with the aim of deriv- 
ing a simple analytical formula for that baseline self-similar 
entropy profile. The hydrodynamical computations in the 
simulations employ entirely different numerical algorithms: 
some use smoothed-particle hydrodynamics (SPH) while the 



others use adaptive-mesh refinement (AMR). Thus, to the 
extent that these two techniques produce convergent re- 
sults, our comparison provides a reliable baseline profile with 
which to interpret the observations. Our comparison shows 
that the codes agree and are presumably reliable outside 
the cores of clusters, where entropy levels are relatively high, 
but disagree inside the cores, where entropy levels are lower. 
This discrepancy is not yet understood and may result from 
differences in how the codes treat small-scale shocks and 
mixing processes. 

While the cosmological parameters used to specify the 
initial conditions in the simulations are similar — ACDM, 
with matter density Qm ~ 0.3, dark-energy density Qa ~ 
0.7, baryon density Q,h ~ 0.04, Hubble constant 2 h ~ 0.7 and 
power-spectrum normalization as ~ 0.9 — the initial condi- 
tions themselves are not identical. The simulations do not 
model the evolution of the same field of density perturba- 
tions. Instead, we compare results for many different clusters 
within representative but not identical volumes of the uni- 
verse, as modeled by each code. 

The paper proceeds as follows. Section 2 establishes a 
non-dimensional framework for comparing the entropy pro- 
files of clusters having different masses then describes the 
two codes we use and presents results from each code. In 
each case we determine the median dimensionless entropy 
profile and the scatter about the median. Knowing the vari- 
ance in entropy owing to the non-steady nature of merging 
should help observers establish whether additional variance 
generated by the stochasticity of non-gravitational processes 
is needed to explain the observations. Section 3 compares the 
results from the two codes and provides a simple analytical 
form for the baseline entropy profile outside the cluster core 
that adequately represents the results of both simulations. 
Because the entropy profiles determined with the two codes 
disagree somewhat in the cluster core, the formula is valid 
only for 0.2 ~ r/r2oo ~ 1-0. Section 4 summarizes our re- 
sults. 



2 SIMULATIONS OF NON-RADIATIVE 
CLUSTERS 

Clusters that form without radiating away any of their ther- 
mal energy are unphysical, but they constitute a useful 
baseline against which to measure the eff ects of radiative 
cooling and non-g ravitational heating (e.g.. lVoit et aljl2002l : 
iFrenk et al J 1 1999ft . In spite of the obvious fact that accre- 
tion shocks generate an enormous amount of entropy dur- 
ing cluster formation. Such clusters are sometimes called 
"adiabatic" clusters in the literature on this subject. Be- 
cause we prefer to reserve the term "adiabatic" for isen- 
tropic processes that do not involve shock heating, we will 
refer to these clusters as "non-radiative" clusters (as in 
iMuanwong et al.ll200 j) . 

This section examines two populations of non-radiative 
clusters generated using both SPH and AMR techniques. 
The properties of these simulated clusters are most easily 
compared if we scale away the dependence of dimensional 
quantities on halo mass, a procedure that would lead to 



1 The classical thermodynamic entropy per particle is s = 
In if 3 / 2 + const, in an ideal monoatomic gas. 



2 We define h = H /(100 krns" 1 Mpc" 1 ). 
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identical temperature, density, and entropy profiles for each 
cluster if clusters were precisely self-similar. Thus, we begin 
by outlining the scaling behavior expected in the self-similar 
case and then remove all the mass- dependent factors when 
analyzing the clusters from each simulation. These simulated 
clusters are not precisely self-similar, in that we typically 
find a ~20% scatter in entropy at a given scale radius outside 
0.1r2oo and a somewhat larger variance inside this radius. 
However, we find no systematic trends in the scaled entropy 
with cluster mass; in other words, we find K(r/r2oo) oc Tx 
to within ~20% across the entire mass range of adequately 
resolved cl usters. Thus, our mod els do not support the sug- 
gestion of IValaeeas et alJ i2003t) that shock heating asso- 
ciated with large-scale structure formation alters the Lx- 
Tx relation by preferentially elevating the entropy levels in 
groups. All of the gravitationally-driven entropy production 
that happens on large scales is well resolved in the simula- 
tions, and the simulated clusters turn out to be self-similar 
over almost two orders of magnitude in mass. 

2.1 Entropy Scaling Laws 

The temperature of a self-similar cluster depends primarily 
on the mass M200 within the radius V2oo, motivating us to 
define the characteristic temperature scale 

GM 20Q ^m p 



T200 



2r2oo 



(1) 



Here and throughout the paper we write the temperature 
in energy units, implicitly absorbing Boltzmann's constant 
into T because observed X-ray temperatures are so often 
quoted in units of keV. Simulations of non-radiative clu sters 
gener ally find that Tx ~ T200 with a scatter of ~10% llVoid 
2005) that apparently depends on the effective resolution of 
a simulation and the numerical techniques it employs. 

The natural entropy scale in the ICM of a non-radiative 
cluster is therefore A200 = T , 200fJe~ 2/ ' 3 , where n e equals 
200f2^ [ 1 times the mean electron density of the universe, 
which would be the mean electron density inside V2oo if the 
electron to dark- matter ratio remained constant. In a ACDM 
cosmology with a baryon density fib = 0.022ft -2 one finds 



n e = (1.45 x 10- 4 cm- 3 )(fi M /0.3)(l + z) 3 , 
Tx ( J21 



giving 



K 



362keVcm z 



IkeV 



/ 1 200 \ 
V Tx ) 



\H{z) 


-4/3 








V0.3 J 



-4/3 



(2) 



where Qm is the current value of the matter-density param- 
eter. Writing A200 in this way makes explicit the fact that 
the observed temperature of a cluster is not necessarily a 
reliable guide to the underlying value of A^oo- If the ICM 
of a real cluster is either hotter or cooler than T200, the 
characteristic temperature of its halo, then one must apply 
the correction factor T200/TX when computing the value of 

A200- 

Radiative cooling introduces another entropy scale into 
the ICM that does not enter into the comparison of simu- 
lated clusters presented here but may be tied to the observed 
K oc T^ 3 scaling in real clusters. Gas of temperature T 
emitting pure thermal bremsstrahlung radiation radiates an 
energy equivalent to its thermal energy in a time period t if 
its specific entropy is 




Figure 1. Dimcnsionless entropy K/K200 as a function of scale 
radius r/r200 for 40 clusters simulated with the SPH code GAD- 
GET. Black squares show the median profile, and the dashed 
line illustrates the power-law relation K/K200 = 1.32(r/r20o) 1 ' 1 • 
Most of the entropy profiles shown lie close to this relation in the 
radial range 0.1 ~ f/»*200 ~ 1-0. At smaller radii, the entropy pro- 
files generally flatten, and their dispersion increases. The shaded 
box shows the range of radii over which the gravity begins to 
depart from a precise inverse square law because of gravitational 
softening. Even though the point at which the entropy profiles 
begin to flatten coincides with the outer edge of this box, we 
suspect that the flattening is real because the better-resolved, 
higher-mass clusters show the same amount of flattening as the 
lower-mass clusters when scaled relative to r200- 



K c « 81keVcirT 



IkeV 



2/3 



t 



14Gyr 



2/3 



(3) 



Because this entropy threshold for cooling is quite similar 

to the observed core entropies of many clusters, it seems a 

2/3 

quite natural explanatio n for the Ap.i oc T x scaling found 
in cluster cores <Voit fc Ponmanll2003ft . However, it is less 
clear why K oc T^ 3 should hold when X»Jf c . 



2.2 SPH Simulations 

The first set of simulated non-radiative clusters we will 
consider was produced by the entropy-conserving ver- 
sion of the smoothed-particle hydrodynamics code G AD- 
GET iSnringel et alJl20oitlSnrin g el & Hernnuistteooll with 
n M = 0.3, n A = 0.7, Q b = 0.045, h = 0.7, and cr 8 = 0.9. 
Most of this set comes from the non-radiative simulation 
described in iKavl {2004) , from which we take the thirty 
most massive clusters, ranging from 2.1 x 10 13 h' 1 M G to 
7.5 x 10 14 ft" 1 M Q . These clusters were modeled within a box 
of comoving length 60 h^ 1 Mpc using an equivalent Plummer 
softening length of 20/i _1 kpc in comoving coordinates. To 
supplement the high-mass end of th e set, we also inclu de the 
ten high-mass clusters simulated bv lKav et alJ i2004l) . using 
the same code and the same cosmological parameters. These 
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Figure 2. Dimensionless entropy K/K200 as a function of halo 
mass M200 at the scale radii r/r2oo = 0.1, 0.5, and 1.0 in 
non-radiative clusters simulated with the SPH code GADGET. 
Squares show entropy at r200i triangles show entropy at 0.5r2oo> 
and diamonds show entropy at 0.1r200- Solid lines give the median 
values of K/K200 at each radius, and dashed lines give the corre- 
sponding values from the power-law relation shown in Figure 111 
No systematic trends with mass are evident. 



clusters, which were extracted from a much larger simula- 
tion volume and then individually resimulated at a spatial 
resolution comparable to that of the 60 h~ x Mq simulation, 
range from 5.5 x 10 14 h' 1 Mq to 8.4 x 10 14 h' 1 Mq. Thus, 
all forty clusters in this overall sample are reasonably well- 
resolved. 

Figure shows the dimensionless entropy profiles of 
these clusters, averaged over radial bins. The average en- 
tropy in each spherical shell is defined to be the mean tem- 
perature in that shell divided by the two-thirds power of 
the mean electron density within the shell. Most of the 
spherically-averaged profiles are virtually identical at r > 
0.1r2oo, consistent with the expectation of self-similarity. 
Within the range 0.2 ~ r/r 20 o ~ 1.0, the power law 
K(r) oc r shown by the dashed line is a good approxi- 
mation, in agreement w ith t he spherical - accre tion models of 
iTozzi fc Normanl fcOO]]) andlVoit et alJ i2003l) and the sim- 
ulat ions~of lBorgarn^et^il] i2002T) . Inside of 0.1r2oo there is 
more diversity. Some of the simulated clusters have nearly 
isentropic cores, while others do not. The flattening of the 
entropy profiles within the core is likely to be a real effect 
because the degree of flattening does not depend on clus- 
ter mass. If the flattening were due to a resolution effect, 
then it would be more pronounced in smaller, lower-mass 
clusters, whose physical size is smaller relative to the reso- 
lution length of the simulation. However, the same kind of 
flattening is s een in the better- resolved, higher-mass clusters 
simulated bv lKav et all (l2004h . 

Figure [5] shows that dimensionless entropy measured at 



a fixed scale radius does not depend on halo mass. No sig- 
nificant deviations from the approximate self-similar profile 
K(r) — 1.32^200 (r/^oo) 1 ' 1 are seen over the entire mass 
range, from 2 x 10 13 h' 1 Mq to 8 x 10 14 h' 1 Mq. We have 
assessed the scatter in dimensionless entropy with the quan- 
tity 

Ait _ K 90 % — K 10 % 
K ~ 



2K 



(4) 



where K x % is the Xth. percentile of the dimensionless en- 
tropy at a given radius. The average value of AK/K in the 
range 0.2 < r/r 200 < 1.0 is 0.12, and over the larger range 
0.02 < r/r 200 < 1.8 its value is 0.28. 



2.3 AMR Simulations 

The other set of clusters we will consider was produced 
by the adaptive-mesh refinement code ENZO fervanlll999l : 
iNorman fc Brvanlll999l : lO'Shea et al.ll2004l) with Q M = 0.3, 
n A = 0.7, n b = 0.04, h = 0.67, and a 8 = 0.9. Twenty- 
one of these non-radiative clusters were modeled within a 
box of comoving length 50/t -1 Mpc using mesh refinement 
to produce an effective resolution of 20/i _1 kpc in comov- 
ing coordinates . Further det ails about this simulation are 
given in lBrvan fc Voit) (120011) . To populate the upper end of 
the mass range, we added ten more clu sters to our s ample , 
drawn from the simulation described in lLoken et alJ J2002fl . 
with a box size of 256 h -1 Mpc and an effective resolution 
of 15 h~ l kpc. That larger simulation assumed a slightly dif- 
ferent cosmology, with Qm ~ 0.3, Qa = 0.7, fib = 0.026, 
h — 0.7, and <jg = 0.928. However, dividing by the appro- 
priate value of K200 when constructing the dimensionless 
entropy profiles compensates for the differing baryon den- 
sity scale. The overall AMR cluster set ranges in mass from 
2.7 x 10 13 /i" 1 M @ to 1.4 x 10 15 h' 1 M©. 

As can be seen in Figure El the entropy profiles mea- 
sured in this sample of simulated clusters are also nearly 
self-similar. Beyond the scale radius 0.2r2oo, the power- law 
slope and normalization of the dimensionless entropy pro- 
files is quite close to those found in the SPH clusters, as 
illlustrated by the dashed line. However, the AMR profiles 
within this radius are significantly flatter than the SPH pro- 
files (§IO- 

Likewise, the dimensionless entropy profiles of our AMR 
clusters also show no systematic trends with mass. Figure 31 
samples the dimensionless entropy in this set of clusters at 
the same three scale radii as in Figure H A gain, the points 
cluster around the median profile, with no dependence on 
M2oo- The scatter in entropy likewise has properties similar 
to those found in the SPH simulation: AK/K = 0.21 over 
the range 0.2 ~ r/raoo ~ 1.0 and AK/K = 0.31 over the 
complete range. 



3 CROSS-COMPARISON OF SIMULATIONS 

It should already be clear from the preceding section that the 
median dimensionless entropy profiles found in both cluster 
simulations are nearly the same outside of 0.2r2oo- Here we 
discuss some of the systematic differences between the two 
sets of simulated clusters and then provide a simple analyti- 
cal formula for the median entropy profile that is consistent 
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Figure 3. Dimensionless entropy K/K200 as a function of scale 
radius r/r200 for 31 clusters simulated with the AMR code ENZO. 
As in Figure most of the entropy profiles shown lie close to 

the relation K/K200 = 1.32(r l /r2oo) 1 ' 1 m t ne radial range 0.2 ~ 
r /^200 ~ 1-0. However, the flattening at smaller radii is more 
pronounced than in the SPH simulation, leading to substantially 
higher entropy levels near the origin. 
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Figure 4. Dimensionless entropy K/K200 as a function of halo 
mass M200 at the scale radii r/r2oo = 0.1, 0.5, and 1.0 in non- 
radiative clusters simulated with the AMR code ENZO. Squares 
show entropy at r200i triangles show entropy at 0.5r200i and dia- 
monds show entropy at 0.1r200- Solid lines give the median values 
of K/K200 at each radius, and dashed lines give the correspond- 
ing values from the power-law relation shown in Figure As in 
Figure |2] there are no discernible systematic trends with mass. 



with both samples. To help observers use this formula, we 
conclude this section with a brief comparison to actual data, 
including a discussion of how offsets of Tx with respect to 
T200 affect the determination of the baseline profile for an 
observed cluster. 



3.1 Entropy within the Core 

The median dimensionless entropy profiles of clusters in our 
two simulation sets agree well outside the cluster core but 
disagree within the cluster core. Figure illustrates the dis- 
crepancy. This type of discrepancy between SPH and AMR 
is nothing new. It was pr eviously hinted at i n the Santa Bar- 
bara cluster comparison (|Frenk et alJll99flh . but that result 
was not definitive because the sample size was a single clus- 
ter. Here we confirm it for a large sample of clusters with 
a range of masses and simulated with substantially higher 
spatial and mass resolution. Thus, we do not attempt to fit 
an analytical form to our median profiles within 0.1r2oo; the 
numerical techniques used to model this region do not yet 
give a reliable answer. 

Because the primary purpose of this paper is to provide 
a baseline profile for observers to use outside the cluster 
core, we leave a detailed analysis of the reasons for this en- 
tropy discrepancy for future work. It is an important prob- 
lem to pursue because of its implications for cooling and 
condensation of gas within cluster cores. Larger amounts of 
entropy production within the core, as in the AMR code, 
will more effectively inhibit cooling there, perhaps mitigat- 
ing the "cooling-flow problem" in clusters of galaxies (see 
iDonahue fc Voit)l2004 for a recent review). 

3.2 Power-Law Approximations 

Previous theoretical work has shown that the entropy 
profiles of non-radiative clusters approximately follow a 
power law with K ( r) oc r 11 (e.g., iTozzi fc Normanl l200ll ; 
iBoreani et"aHl200dlVoit et alJl2003l) . but these efforts have 
not provided a normalization for that power-law profile in 
a form that is useful to observers. Here we rectify that sit- 
uation. If we fix the power-law slope of the entropy pro- 
file at 1.1 and fit the SPH clusters in the radial range 
0.2 < r/r2oo ^ 1.0, we find 



K(r) = 1.32 ± 0.03 /Goo (r/r 20 o) 11 ■ 
Doing the same for the AMR clusters yields 
K{r) = 1.41 ± 0.03 #200 (r/raoo) 1 ' 1 , 



(•>) 



(6) 



a normalization just slightly higher than that for the SPH 
clusters. 

Figure |S] compares the power-law fit from equation 
with the median profiles from both simulations and with two 
other analytical entropy profiles. One of the analytical pro- 
files is constructed from t he analytical temp erature and den- 
sity profiles developed bv lRasia et alJ i2004) to fit their SPH 
models of non-radiative clusters. The other is the entropy 
profile of intracluster gas with an NFW density distribution 
with concentration c = 5 assuming that it is in hydrostatic 
equilibrium within a dark-m atter halo whose density distri- 
butio n has the same shape (IVoit fc BrvanllioOll : IVoit et ail 
2002). All of these profiles agree very well with the power- 
law profile of equation (jy in the range 0.2 ^ r/r^oo ^ 1-0. 
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Figure 5. Median entropy profiles from cluster simulations with- 
out cooling or non-gravitational heating. A solid line connecting 
crosses shows the median profile for the SPH clusters. A dashed 
line connecting squares shows the median profile for the AMR 
clusters. Error bars give the 10% percentile to 90% percentile 
range. The long dashed line illustrates the power-law approxima- 
tion K(r) = 1.32A"200 (?"A"20o) ■ The dotted line gives an an- 
alytic al entropy profile derived from simulations by iRasia et alj 
(2004). The dot-dashed line shows the profile i^NFwM corre- 
sponding to an NFW gas-density profile of concentration c = 5 
that is in hydrostatic equilibriu m within a dark-matter density 
profile of identical shape iVoit et al.ll2002l) . Notice that both of 
the median profiles agree very well in the range 0.2 < r/r 2 oo < 1-0 
and that the analytical approximations accurately represent the 
median profiles in this range. However, the AMR and SPH me- 
dian profiles differ by as much as a factor of two within the cluster 
core. 

Also, the power-law profile remains a good representation 
of our SPH simulat ions down to ~ 0.1r2oo. We note, how- 
ever, that the IRasia et al.l <l2004t) profile shows much less 
flattening in the core than the simulations analyzed in this 
paper, perhaps because it is based on an earlier version of 
GADGET that does n ot explicitly conserve entropy (see also 
lAscasibar et al.ll2003l) . 

Up to this point, we have been fitting the median pro- 
files with a K(r) oc r 1 ' 1 power law because that is the stan- 
dard power-law index in the literature, and it appears to be 
consistent with the hig hest quality cluster observations (e.g., 
iPratt fc Arnaudll2003l) . However, Figure © suggests that a 
power-law index of 1.1 might be slightly too shallow to be the 
best representation of the median profiles outside of 0.2r2oo- 
Fitting a power-law of index 1.2 to the median profiles gives 
the following results: 

K(r) = 1.43 ±0.01 ivaoo (r/raoo) 1 - 2 . (7) 
for the SPH clusters and 

K(r) = 1.48 ± 0.02 K 200 (r/r 2m f- 2 (8) 

for the AMR clusters. Figure |7] shows that dividing the me- 
dian profiles by equation J7J makes the resulting profiles 
flatter, indicating that K{r) oc r 1 ' 2 is a better description 



Figure 6. Median dimensionless entropy profiles from 
simulations, divided by the power-law profile K(r) = 
1.32 A^oo (x/f20o) ■ Crosses connected by a solid line show the 
median profile from the SPH simulation, and squares connected 
by a dashed line show the median profile from the AMR simu- 
lation. Dashed and dot-dashed lines show the IRasia et al.l i2004h 
and NFW-like profiles, respectively. A slight rise in the median 
points as radius increases beyond 0.2r2oo suggests that the as- 
sumed power-law index of 1.1 is slightly too small. 

of the outer parts of non-radiative clusters. Indeed, if we al- 
low the power-law index to be a free parameter and fit the 
median profiles in the range 0.2 ^ r/r 2 oo ^ 1-0, we find 

K(r) = 1.45 ±0.01 K 200 {r /r 2m f- 21±om . (9) 

for the SPH clusters and 

K(r) = 1.51 ± 0.03 K 200 (r/r 200 ) 1 - 24±om (10) 

for the AMR clusters. In all of these fits the error bars cor- 
respond to la. 

3.3 Applications of the Baseline Profiles 

The entropy profiles computed here for non-radiative clus- 
ters provide a baseline for assessing the impact of non- 
gravitational processes on the intracluster medium. How- 
ever, in order to compare the entropy profiles of a real cluster 
to these self-similar baselines, one needs to know the value 
of -K200 for the cluster. This characteristic entropy scale can 
be simply computed from equations (1) and (2) if the clus- 
ter mass M200 has been accurately measured. Otherwise, 
one must infer K200 from Tx and a relation between T200 
and Tx or, equivalently, a relation between M200 and Tx. 
Figure [S] compares entropy values measured at 0.1r2oo with 
those predicted for self-similar clusters. The long-dashed and 
dotted lines show the predictions for non-radiative AMR and 
SPH clusters, respectively. Here we set Tx = T200 because 
any offset between Tx and T200 is small compared with the 
differ ence between t he t wo s imulation se t s. As shown by 
iPonman et"aD (120031) and lVoit fc Ponmanl ll2003l) . the mea- 
surements clearly do not agree with the self-similar models, 
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Figure 7. Median dimensionless entropy profiles from 
simulations, divided by the power-law profile K(r) = 
1.43^200 { r l r 20o) 1 '' 2 ■ Crosses connected by a solid line show the 
median profile from the SPH simulation, and squares connected 
by a dashed line show the median profile fro m the AMR simula - 
tion. Dashed and dot-dashed lines show the iRasia e t al. (2004) 
and NFW-like profiles, respectively. A power-law index of 1.2 
seems to be a better description of the entropy profile beyond 
0.2r200 than the standard index of 1.1. 

which predict that K{0.1r2oo) ex Tx. Instead, the measure- 
ments track the cooling threshold K C (T) oc T^/ 3 . We wish to 
point out, however, that the hot clusters are consistent with 
non-radiative m odels at this rad i us. Th is finding contrasts 
with Figure 1 of lVoit fc Ponmanl J2003), in which the locus 
for self-similar clusters is mistakenly a factor of two too low, 
owing to a unit conversion error. 

Figure [5] compares entropy values measured at 7-500 ~ 
0.66r2oo, within which the mean mass density is 500 times 
the critical density, with the baseline profiles. The solid line 
shows the baseline entropy level derived assuming Tx = 
T200, which slightly exceeds the measured entropy levels at 
this radius in hot clusters. This apparent shortfall in the 
observed entropy levels of hot clusters goes away when we 
account for the difference between Tx and T2oo- A dashed 
line shows the baseline entropy level at rsoo computed us- 
ing values of .K 2 00 de rived from the M200-TX relation of 
ISanderson et alJ J2003f) . and this line is consistent with the 
entropy measurements at rsoo in hot clust ers. This consis- 
tency again contrasts with the results of IVoit fc Ponmanl 
(2003), in which the locus for self-similar clusters was mis- 
takenly placed too low. However, clusters below about 6 keV 
still show a clear entropy excess, which is even more pro- 
nounced when the observed M200-TX relation is used to com- 
pute 7\200- 



4 SUMMARY 

Our intention in this paper has been to provide a simple ana- 
lytical form for the entropy profiles of non-radiative clusters 




T x (keV) 

Figure 8. Relationship between core entropy and the cooling 
threshold. Each point with error bars shows the mean core en- 
tropy -Ko.li measured at 0.1r2oo, for eight clusters within a given 
temperature bi n, and small crosses s how measurements for indi- 
vidual clusters (Ponman et al. 2003). The long-dashed and dot- 
ted lines show self-similar relations calibrated using the median 
values of K0.1 derived from our AMR and SPH simulations, re- 
spectively, assuming that Tx = T2oo- High-temperature clusters 
appear to be consistent with the median self-similar profiles, but 
the trend to lower-temperatures more closely tracks the cooling 
threshold K C (T) (solid line), defined to be the entro py at which 
the cooling time equals 14 Gyr <Voit fc Ponmanl2003l) . The short- 
d ashed line shows the predicted entropy at 0.1r2oo in th e model 
of lVoit fc BrvarJ teOOllV 

to use as a baseline when trying to measure the impact of 
non-gravitational processes on the intracluster medium. To 
that end, we analyzed two different sets of simulated clus- 
ters, one created with a Lagrangian SPH code and the other 
with an Eulerian AMR code. Thirty to forty entropy profiles 
were produced by each code and these profiles were found to 
be approximately self-similar, with the K(r/r2oo) oc Tx scal- 
ing expected of non-radiative clusters. The simulated profiles 
depend very little on halo mass once the expected scaling 
is divided out. This result confirms that non-gravitational 
processes are necessary to produce the observed scaling re- 
lations of clusters. 

The median entropy profiles from the two simulations 
agree to within 7% outside of 0.2r2oo but disagree in the 
cluster core. In the outer parts of clusters the power-law 
profile K(r) = 1.32 .ff 200 (r/raoo) 1 ' 1 is a good representation 
of the baseline profile expected in the non-radiative case. 
However, our results suggest that the baseline profile in the 
radial range 0.2 ^ r/r2oo ^ 1.0 is better fit by a K(r) tx r 1 ' 2 
po wer law, rather than th e standard K (r) oc r 1 ' 1 law found 
bv lTozzi fc Normaiil (1200 ll) . Inside of 0.2r2oo the discrepancy 
between the AMR clusters and the SPH clusters is quite 
substantial. Tracking down the origin of this discrepancy is 
important, because radiative cooling rates in clusters depend 
on the core entropy level. 

Our comparison between the self-similar entropy pro- 
files derived from these simulated clusters with measure- 
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Figure 9. Entropy at rsoo as a function of cluster temperature. 
Each point with error bars shows the mean value of A" (rsoo) 
implied by the density and tem perature profiles of eight clus- 
ters within that temperature bin I Ponman ct al . 2003), and small 
crosses show measurements for individual clusters. The solid line 
shows the median entropy at rsoo f° r self-similar clusters, assum- 
ing that Tx = T2oo- The dashed line shows how the mapping of 
this median entropy onto Tx changes when the Tx(M2oo) rela- 
tionship observed bv lSanderson et alJ 120031) is used to determine 
Tx. The most massive clusters are consistent with the self-similar 
clusters modeled without non-gravitational processes when this 
observational M200-TX relation is used. 



merits of entropy at 0.1 r2oo and rsoo in real cluster s updates 
and corrects the findings of IVoit fc Ponman! (2003). There is 
a clear entropy excess in cool clusters, presumably stemming 
from non-gravitational processes. However, clusters hotter 
than ~ 6keV appear to converge to the self-similar profile 
at radii ~ 0.1 r2oo- 

We thank Trevor Ponman for persistently asking the 
questions that motivated this project and and Volker 
Springel for generously allowing us to use GADGET2 before 
its public release. The SPH simulation data were generated 
using COSMA, the 670-processor COSmology MAchine at 
the Institute for Computational Cosmology in Durham, as 
part of the programme of the Virgo Supercomputing Consor- 
tium. GMV received support for this project from NASA's 
Astrophysics Theory program through grant NNG04GI89G. 
STK is supported by the UK Particle Physics and Astron- 
omy Research Council (PPARC). GLB thanks PPARC and 
the Leverhulme Trust for support. 



REFERENCES 

Allen, S. W., & Fabian, A. C. 1998, MNRAS, 297, L57 
Arnaud, M., & Evrard, A. E. 1999, MNRAS, 305, 631 
Ascasibar, Y., Yepes, G., Muller, V., & Gottlober, S. 2003, 

MNRAS, 346, 731 
Balogh, M. L., Pearce, F. R., Bower, R. G., & Kay, S. T. 

2001, MNRAS, 326, 1228 



Borgani, S., Governato, F., Wadsley, J., Menci, N., Tozzi, 
P., Quinn, T., Stadel, J., & Lake, G. 2002, MNRAS, 336, 
409 

Borgani, S., Finoguenov, A., Kay, S.T., Ponman, T.J., 
Sprin gel, V., Tozz i, P., & Voit, G.M. 2005, MNRAS, in 
press, astro-ph/0504265 

Bryan, G. L. 1999, Computing in Science and Engineering, 
1, 46 

Bryan, G. L., & Voit, G. M. 2001, ApJ, 556, 590 

David, L. P., Jones, C, & Forman, W. 1995, ApJ, 445, 578 

Donahue, M., & Voit, G. M. 2004, in Carnegie Astrophysics 

Series, Vol. 3: Clusters of Galaxies: Probes of Cosmologi- 

cal Structure and Galaxy Evolution, ed. J. S. Mulchaey, 

A. Dressier, & A. Oemler, in press 
Edge, A. C, & Stewart, G. C. 1991, MNRAS, 252, 414 
Eke, V. R., Navarro, J. F., & Frenk, C. S. 1998, ApJ, 503, 

569 

Evrard, A. E., & Henry, J. P. 1991, ApJ, 383, 95 

Frenk, C. S., et al. 1999, ApJ, 525, 554 

Kaiser, N. 1986, MNRAS, 222, 323 

— . 1991, ApJ, 383, 104 

Kay, S. T. 2004, MNRAS, 347, L13 

Kay, S. T., Thomas, P. A., Jenkins, A., & Pearce, F. R. 

2004, MNRAS, 355, 1091 
Kay, S. T., Thomas, P. A., & Theuns, T. 2003, MNRAS, 

343, 608 

Loken, C, Norman, M. L., Nelson, E., Burns, J., Bryan, 

G. L., & Motl, P. 2002, ApJ, 579, 571 
Markevitch, M. 1998, ApJ, 504, 27 

Muanwong, O., Thomas, P. A., Kay, S. T., Pearce, F. R., 
& Couchman, H. M. P. 2001, ApJ, 552, L27 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1995, MN- 
RAS, 275, 720 

Norman, M. L., & Bryan, G. L. 1999, in ASSL Vol. 240: 

Numerical Astrophysics, 19 — h 
O'Shea, B. W., et al. 2004, |astro-ph/0403044| 
Ponman, T. J., Sanderson, A. J. R., & Finoguenov, A. 2003, 

MNRAS, 343, 331 
Pratt, G. W., & Arnaud, M. 2003, A&A, 408, 1 
Pratt, G. W., & Arnaud, M. 2005, A&A, 429, 791 
Rasia, E., Tormen, G., & Moscardini, L. 2004, MNRAS, 

351, 237 

Sanderson, A. J. R., Ponman, T. J., Finoguenov, A., Lloyd- 

Davies, E. J., & Markevitch, M. 2003, MNRAS, 340, 989 
Springel, V., & Hernquist, L. 2002, MNRAS, 333, 649 
Springel, V., Yoshida, N., & White, S. D. M. 2001, New 

Astronomy, 6, 79 
Tornatore, L., Borgani, S., Springel, V., Matteucci, F., 

Menci, N., & Murante, G. 2003, MNRAS, 342, 1025 
Tozzi, P., & Norman, C. 2001, ApJ, 546, 63 
Valageas, P., Schaeffer, R., & Silk, J. 2003, MNRAS, 344, 

53 

Valdarnini, R. 2003, MNRAS, 339, 1117 

Voit, G. M. 2005, Rev. Mod. Phys., 77, 207 

Voit, G. M., Balogh, M. L., Bower, R. G., Lacey, C. G., & 

Bryan, G. L. 2003, ApJ, 593, 272 
Voit, G. M., & Bryan, G. L. 2001, Nature, 414, 425 
Voit, G. M., Bryan, G. L., Balogh, M. L., & Bower, R. G. 

2002, ApJ, 576, 601 
Voit, G. M., & Ponman, T. J. 2003, ApJ, 594, L75 
White, S. D. M., & Frenk, C. S. 1991, ApJ, 379, 52 



